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Abstract 



We consider a ring of strongly interacting electrons connected to two ex- 
ternal leads by tunnel junctions. By studying the positions of conductance 
resonances as a function of gate voltage and magnetic flux the interaction pa- 
rameter g can be determined experimentally. For a finite ring the minimum 
conductance is strongly influenced by device geometry and electron-electron 
interactions. In particular, if the tunnel junctions are close to one another 
the interaction-related orthogonality catastrophe is suppressed and the valley 
current is unexpectedly large. 
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Typeset using REVT^ 



When the size of an electronic system is reduced, a rich variety of new "mesoscopic" 
phenomena becomes experimentally observable. Some of the new phenomena are essentially 
classical, owing their existence to the granularity of the electric charge and the system 
size -dependence of various energy scales. The most widely studied example of this type 
is Coulomb blockade. A different category of mesoscopic phenomena is entirely quantum 
mechanical, and is due to the fact that the phase coherence length at low temperatures is 
comparable to the system size, giving rise to a number of interference effects. A particular 
example is the existence of a persistent current in the ground state of a mesoscopic ring. In 
this Letter we study the interplay between two mesoscopic phenomena. Coulomb blockade 
and Aharonov-Bohm interference, using an exactly solvable model. 

We consider a system consisting of a small ring of interacting electrons connected to 
two non-interacting leads by tunnel junctions. The tunnel junctions are at positions Xl and 
xr, respectively. The ring is capacitively coupled to an external gate electrode and may be 
pierced by a magnetic flux. We consider a small AC voltage applied to the right lead and 
wish to evaluate the current at the left junction. A straightforward apphcation of Kubo 
formula |I| yields the current 

{hit)) = |_'dt'e-^(*'-*)Tr{pG [lR{t'), h{t)]} 

where V{t) = Vqc"*^* is the applied voltage and pc is the equilibrium density matrix. The 
quantity on the right hand side is recognized as the retarded current- current correlation 
function. It is most readily evaluated in imaginary time. 

The Hamiltonians for the leads and the connection between the leads and the ring are 
given by 

HL = Y,^kLc{{k)cL{k) (1) 

k 

HlT = tLc\{xL)i'{xL) + tlij\xL)CL{xL) (2) 

where cl are operators on the left lead and if) are operators on the ring. The Hamiltonians 
for the right lead are defined analogously. Thus, we take the leads to be non-interacting 
and couple them to the ring with tunnel junctions at positions xl and xr. To calculate 
the current-current correlation function we define the generating functional Z[Ji, Jji\ = 
Tt exp[—/3H — i jj^^ dT{JL{T)lL{T) + Jr{t)Ir{t))] and integrate out the free fermions in the 
leads. That yields 

Z[Jl, Jr] = ZiZHTrexp|-/5ifring - w lo^ ^r' 

- f JL(r))(l + fJL{r')W{T, XL)GL{r - r'; x^, Xi)^(r', xl) (3) 



+ ^{1 + '-^ 



jR{r)){l - f J^(r'))^t(^, XR)GR{r - r'; xr, XR)^(r', xr)\ } 



where Z^jr and Gl/r{t] x, x') are the partition functions and free fermion propagators in the 
leads, and -faring is the Hamiltonian for an isolated ring. The imaginary time ordered correla- 
tion function x(ri — T2) = —{Tt-{Il{ti)Ir{t2))) is obtained by differentiating Z[Jl, Jr\ with 
respect to Jl and Jr. To simplify the notation we introduce the four-operator expectation 
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value A{ti,T,T2,t') = A{ti-T2,Ti~T,T2-t') = {Tr{lp^{T,XL)lp{Tl,XL)lp^{T2,XR)i!{T',XR))) 

and its Fourier transform A{iuJn,iuJn,iuj'^)- To the lowest nonzero order in the tunnehng 
matrix elements x is given by 

1 r 

X [GL{iuJn - i^n, xl,xl) - G L{iuj„; X L, X l)][G R^iu'^; X R, X r) - GR{iLj'^ + iVtn] xr,xr)] 

(4) 

If we regard the leads as infinite (rather than semi-infinite), the propagators in the leads are 
easy to evaluate and yield GLiitOn', xl, xl) = —i^DL{eF)siga.{uJn) where -Di(e) is the density 
of states in the left lead. 

To evaluate the four-operator product A(ia;„, ia;^, ia;^) we must specify the Hamiltonian 
for the ring. We choose to work with the simplest exactly solvable interacting model, the 
spinless Luttinger model. In the bosonized form the Hamiltonian reads 

H,.n, = ^[-{N - Nof + gv{J - JoY] + J2 ^^^1^1^% 

where and J are zero modes associated with the total charge and total current. Since 
the numbers of clockwise and counterclockwise moving electrons on the ring must both be 
integers, the quantum numbers N and J must satisfy (—1)^ = (— I)"'- The gate voltage and 
magnetic flux determine the parameters Nq = GVg/e and Jq = 2$/$o which in turn deter- 
mine the ground state charge and current. The parameter g is a measure of the interaction 
strength, and equals one for non-interacting electrons For future use we also define the 
shorthand notation 7 = \{g + g~^) — 1 which vanishes in the non-interacting limit. 

Due to time-ordering the exact expression for A is quite complicated although in principle 
straightforward. Since we consider only the lowest order in the tunneling Hamiltonian, our 
analysis is valid only sufficiently far from the conductance resonances. Therefore, we can 
assume that the ground state is separated from the excited states by an energy gap 5e 
that is larger than fc^T. This approximation is basically similar to the one used by Fazio 
and co-workers |Q for an interacting ring connected to superconducting leads. We also 
neglect events with all four imaginary times approximately equal since their contribution is 
negligible. That allows us to evaluate A in an approximate fashion, and gives 

A{iVtn, -iujn, -i^^'n) ~ h(3[5n„fiG{iuJn, Q)G{iuj'^, 0) 

-5f7„,a;„-c^;G'(zCJ„, XL - XR)G{iuJn - i^n, Xr - Xl)] 

where G{t,x) = — (TT-(?/'(r, x)-?/'^(0, 0))) is the Green's function for interacting electrons on 
the ring and G{iuJn,x) is its Fourier transform with respect to the imaginary time variable 
T. For non-interacting electrons we can apply Wick's theorem and find that this expression 
is exact. Substituting this into the expression @ gives, after proper analytic continuations, 
the DC conductance 



a ^ j^-^^DLieF)DR{eF) T diu (-^^ ) G-\u;,xl - xr)G-'^{u,xr - xl) 



00 



^ J~^^L^'F)DR{eFW''\u = 0,xl-xr)\' 



(5) 
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where the last expression is vahd if the system is further than fc^T from a resonance. 

This expression can also be understood using a scattering matrix approach regarding the 
ring as a (complicated) scatterer for the free electrons in the leads as discussed in a specific 
case by Jagla and Balseiro . From that point of view our basic approximation is that one 
scattering event is completed before another one takes place — the approximation breaks 
down near resonance when the dwell time for extra electrons in the ring is large. 

Now we turn to evaluating the retarded Green's function for interacting electrons in 
a finite ring at a finite temperature. We use the low-energy expansion iIj{x) = ip-\-{x) + 
tp^{x) where ip±{x) are clockwise and counterclockwise moving fermions so that G{t,x) = 

G++{t,x) + G (r, x). Following Haldane 0, the correlation functions can be evaluated 

exactly and we obtain (p = ±1) 



l?3(^g-laAro-x^^,e-^g"'")l?3(^gapJo-XJ,e-^g")+1?2(^g~^aA^O-a:iv,e-23~'")l?2fa^^ 

i93(ig-iaAfo,e-29-ia)^3(igQ,pJg,e-29")+i92(jg~i«Afo,e-29~"'")i?2(j9Q:p^0.e~^''") ^ ' 

l^'i(0,e-")r+^ 
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where 7 = }^{g + g^^) — 1 and we introduced the shorthand notation a = ZE!^^ = 
j^{ig~^vT — px), and xj = j^{igvT — px) . Here a is a short distance cutoff for the interaction, 
and is of the order of the lattice spacing. The Jacobi theta functions -i^s and '(?2 arise from 
the g = modes with N and J both even or odd, respectively, and the -^i arises from the 
bosons with q 0. The appearance of doubly periodic elliptic functions is natural since the 
Green's function must be periodic in x and antiperiodic in r. The Jacobi functions appear 
also in the partition function of an isolated ring (essentially the x- and r-independent parts 
of the above expression) and therefore in most equilibrium properties of mesoscopic rings 
like persistent currents 

The parameters that are most readily accessible in an experiment are the gate voltage 
and the magnetic flux. They enter only the q = part of the Hamiltonian which we can 
re-write as Hq = \Ec{N — N^f + ^^{J — JqY where vp = gv is the Fermi velocity of a 
non-interacting system with the same density and = is the charging energy. The 
conductance resonances correspond to values of the gate voltage and magnetic flux at which 
the ground state quantum numbers and J change (degenerate ground state). Hence, 
the positions of conductance resonances can be determined from a simple charging energy 
model with a single-particle Aharonov-Bohm term — note, however, that is not simply 
given by the geometric capacitance (it is non-zero even for non-interacting electrons). In 
the (V^, <l')-plane the resonance positions form a network the shape of which depends on 
the interaction parameter g. We suggest therefore that the interaction parameter can be 
experimentally measured by studying the trajectories of conductance maxima as a function 
of the gate voltage and magnetic flux. For non-interacting systems the resonance positions 
form a lattice of diamond shaped parallelograms, whereas for repulsive interactions [g < 1) 
there are some values of the gate voltage such that a resonance condition is not met for 
any $ as indicated in Figure 0. For attractive interactions, there are ranges of $ such that 
the total current in the ring, J, is independent of Vg. In that case electrons can tunnel 
into and out of the ring only as pairs of clockwise- and counterclockwise movers, which is 
reminiscent of Cooper pair tunneling through a superconducting grain [Q. From now on we 



4 



only consider repulsive interactions. 

We use the expressions (H) and @ to analyze the dependence of the conductance on the 
external parameters. The dependences on Vg and $ are qualitatively similar and in Figure 
H we show the conductance as a function of Vg. Near a resonance we obtain the limiting be- 
havior |G"''^*(0, x) p ~ C /{def where 5e is the energy cost of changing the number of electrons 
in the ring by one (we still assume 5e > ksT). In a generic case a resonance corresponds to 
a degeneracy for the addition or removal of either a clockwise or counterclockwise moving 
particle, and the prefactor C is independent of device geometry. However, since we have two 
control parameters Vg and $, we can use them to bring both clockwise and counterclockwise 
modes to resonance simultaneously. These special parameter values correspond to slope 
changes in the trajectories of conductance resonances in Figure |l]. Since at these double 
resonances the Green's function has significant contributions from wave vectors kF+ and 
—kp-, its absolute square has components with wave vector kp^ — {—kp-) = 27iN/L. Thus, 
we find that near a double resonance the amplitude of the conductance maximum varies 
periodically with the separation between the tunnel junctions as (1 + cos[^^(xl — x/j)]). 
Away from a double resonance one channel dominates and the amplitude of these oscilla- 
tions is exponentially small. Since the wave vector of these oscillations is essentially 2kF, 
they can be observed only in low-density systems or using local probes like STM. The in- 
terference effects are smeared out by temperature which leads to a different temperature 
dependence of the conductance for different device geometries near a double resonance: if 
N\xl — xrI ~ nL, where n is an integer, the conductance decreases with temperature due 
to reduced interference, whereas for N\xl — xr\ ~ (n + i)L the conductance increases with 
T. 

From the expression (^ we see that the conductance is reduced due to interactions by a 
factor (a/L)^"^ which can be attributed to an orthogonality catastrophe that has previously 
been studied in the context of quantum dots in the fractional quantum Hall regime |jlO[- The 



exponent governing the resonance line shape for small 6e, a ~ l/(5e)^, is independent of the 
interaction parameter g, which is a consequence of a finite minimum energy of the bosonic 
modes. The resonant contribution dominates for 6e ^ 6ec = sin ^ZlliiijMl^ (up to 

logarithmic corrections); for 6e ^ 6ec the valley conductance levels off to a constant value 
proportional to {a / \x l — x rIY'^ as seen in Figure ^ For large separations Ax = \xl — xr\ the 
crossover point 6ec exceeds half of the resonance spacing and the crossover is not observed. 
The two limiting behaviors can be combined to give the approximate line shape 



a{6e) ~ TlTr 



{6ej^ \ l) ^ [2^) 



L ( ttAx 
sm 



V2a V ^ 



g^log|sm(^)| 
7 



2 



(7) 



where 7 = |((7+(7^^) — 1 and F^^/i? = |iL/_Rp-DL//j(ei7') are the line widths for a non-interacting 
system. The last factor gives rise to a logarithmic dependence on Ax in the non-interacting 
limit. 

At T = the smallest 8e that we can consider is determined by when terms that are 
higher order in the tunneling Hamiltonian become significant. Due to the orthogonality 
catastrophe that happens when the first term in (|^) is of order unity. Therefore, we expect 
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that at T = the height of the resonance peak is independent of g, and the peak width is 
given by 6e such that TlTr (\/2-|) (5e)~^ 1. A simple Breit-Wigner formula a ~ (^^p^p2 
gives a peak-to- valley ratio 2(r/A)^ where A is the separation between adjacent resonances 
and r is the resonance width. In the present case this simple connection does not hold: the 
line width at T = is reduced by a factor [a/L)"' whereas the valley current (for small Ax) 
is only suppressed by factor {a/AxY'^. The former suppression factor can be identified as 
the life time of a charged excitation of the ring while the second one is the off-resonance 
probability of transmission through the ring. The valley current is therefore anomalously 
large for small Ax since effectively the system size is replaced by Ax and the orthogonality 
catastrophe is less severe. Since at T = the peak widths are reduced by interactions, the 
resonance peak heights at a finite temperature are suppressed due to thermal broadening. 
Thus, the main effect of interactions is to change the peak-to-valley ratio in a way that 
depends on device geometry and temperature. 

The experimental possibilities for the study of nanostructures like the one we consider 
are developing rapidly. New techniques like conducting organic molecules and carbon com- 
pounds are emerging to complement the conventional semiconductor structures. In par- 
ticular, it was recently demonstrated [|lT| that carbon nanotubes exhibit coherent electron 



transport and can be used to fabricate nanoscale ring structures. We believe these devices 
can be used to experimentally study the system we have analyzed. 

In conclusion, we have considered tunneling through a finite strongly interacting system 
within the framework of an exactly solvable model. We find that the positions of conductance 
resonances in the (V^, $)-plane can be used to determine the interaction parameter g. We 
conclude that at T = the heights of resonance peaks are unaffected by interactions but 
due to the narrowness of T = resonances, the peak conductance at a finite temperature 
is reduced by interactions. The valley current depends on both interactions and the device 
geometry. Near a double resonance we find that the heights of resonance peaks depend on 
device geometry due to interference between different current carrying processes. 
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FIGURES 




FIG. 1. Positions of conductance resonances in the (V^,$)-plane for interacting electrons 
(repulsive interactions, g = l/\/2). The labels N and J denote the ground state charge and 

current as a function of the external parameters, and the shaded area indicates the domain of 
validity of our analysis. The line segments with different slopes correspond to fluctuations in the 
numbers of clockwise and counterclockwise moving electrons, respectively. Inset: device geometry. 
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FIG. 3. Conductance vs. \xl — xr\ near a conductance minimum for interaction parameters 



g = 0.9, g = 0.7, and g = 0.5. The temperature is T = 0.1^. 
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